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Abstract 



o 



Achromatic breaks in afterglow light curves of gamma-ray bursts (GRBs) arise naturally if the product of the jet's Lorentz factor y 
and opening angle ®j satisfies y®j » 1 at the onset of the afterglow phase, i.e., soon after the conclusion of the prompt emission. 
Magnetohydrodynamic (MHD) simulations of collimated GRB jets generally give y&j < 1, suggesting that MHD models may be 
inconsistent with jet breaks. We work within the collapsar paradigm and use axisymmetric relativistic MHD simulations to explore 
the effect of a finite stellar envelope on the structure of the jet. Our idealized models treat the jet-envelope interface as a collimating 
rigid wall, which opens up outside the star to mimic loss of collimation. We find that the onset of deconfinement causes a burst of 
acceleration accompanied by a slight increase in the opening angle. In our fiducial model with a stellar radius equal to 10 4 ' 5 times 
that of the central compact object, the jet achieves an asymptotic Lorentz factor y ~ 500 far outside the star and an asymptotic 
opening angle ; a* 0.04 rad 2°, giving y@j ~ 20. These values are consistent with observations of typical long-duration GRBs, 
and explain the occurrence of jet breaks. We provide approximate analytic solutions that describe the numerical results well. 



20091) . Such thermally-driven jets accelerate be cause of ex 
pansi on of t h e relativistically hot g a s (e.g., Zhang etal 
20031 l2004t iMorsonv et ail 120071: IWangetall 120081) 



However, self-consistent simulations including neutrino 
physics have not yet succeeded in launching relativistic jets 
ilNagataki et ail 120071: iTakiwaki et al l 120091: iNagatakl 120091: 



however, see lHarikae et al.M2009f) . The inclusion of magnetic 



fields enables the extraction of rotational energy from the 
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i 1. Introduction 

o ■ 

■ Relativistic jets are ubiquitous features of many accreting 
j3 black holes and neutron stars. They are found in systems 
l— ~ 1 spanning an enormous range of compact object mass: from x- 
ray binaries (XRBs) and gamma-ray bursts (GRBs) with neu- 
j>. tron stars and stellar-mass black holes (BHs) of mass M ~ 
I |1.4-2OM , to active galactic nuclei (AGN) with BHs of mass 
•M ~ 1O 6 -1O 1O M . The physics of jet production appears to be 
robust and insensitive to the details of the central object. 

Two characteristic features of relativistic jets are the fol- 
lowing: (i) Jets are accelerated efficiently and achieve large 
Lorentz factors y, ranging to above several hundred in the case 
of GRBs. (ii) Jets are strongly confined, with opening angles 
&j < 0.1 rad. It is generally assumed that these two properties 
. £^ are related, though the exact relation is not well understood. 

Long-duration GRBs (which we hereafter refer to simply as 
GRBs) are particularly interesting in this regard since we have 
measurements of both y and ® , ( Piranl 120051) . Many GRBs 
have y > 400 (ILithwick & Saril 1200 lb . with jet opening an- 
gles broadly dist ributed around a typical value ®j ~ 0.05 rad 
(iZeh et all 120061) . These estimates imply that y®j ~ 10-30. 

Jet acceleration requires an energy source. In the ab- 
sence of magnetic fields this source may be thermal energy 
injected into the jet by annihi l ating neutrinos from an ac- 
cretion disk dKohri et all 120051: jChen & Beloborodovl 12007 
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TMNOi). In such mod- 



els, acceleration is caused by the expansion of magnetic 
fields which have been twisted and amplified by th e rota 



tion of the central co mpact obj ect (IBegelman & Lil 11 994 
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acceleration mechanism in the present paper. 
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tic magnetohydrodynamic (MHD) jets confined by an external 
medium should have y® j < 1 . Their result appears to be in 
conflict with the GRB data summarized above. On the other 
hand, a completely unconfined split-mono pole flow achieves 
larger opening angles and gives y®j > 10 JtMN09I) . How do 
these results relate to GRBs? According to the c ollapsar model 
dWooslevl Il993t iMacFadven & Wooslevl Il999t) . GRB jets are 
produced by an accreting stellar-mass black hole at the cen- 
ter of a collapsing massive star. As the jet propagates out, it 
is collimated by the pressure of the shocked stellar envelope. 
However, once the jet emerges from the star it is effectively in 
vacuum, since the pressure of the external interstellar medium 
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or stellar wind is far lower than the internal pressure in the jet. 
Thus a GRB jet corresponds to a hybrid scenario which involves 
initial confinement followed by free propagation. What are the 
asymptotic properties of such a jet far from the star? In particu- 
lar, do es yQj have a value of order unity, as in confined models 
Hoi), or is it com parable to the values seen in unconfined split- 
monopole models jTMN09t) ? 

We have carried out high-resolution simulations that model 
the behavior of a relativistic magnetized jet emerging from stel- 
lar confinement. Our ideal MHD jets are confined out to a cer- 
tain distance by a rigid wall of the appropriate shape to mimic 
the confining effect of the shocked stellar envelope (see 
We then allow the jets to be effectively unconfined beyond 
this distance. The simulations are meant to represent quasi- 
steady GRB jets, i.e., we ima gine that an initial (p ossibly non- 



relativistic) precursor jet (e.g jBurrows et al.ll2007l) has already 



cleared the way for an ultrarelativistic magnetized jet which 
later powers the GRB. It is the latter jet that our simulations 
are meant to represent. The numerical setup is discussed in $2] 
and the results are described in 50 We conclude in ^4] with a 
discussion of the results in the context of GRBs. We work with 
Heaviside-Lorentzian units, neglect gravity, and set the speed 
of light and the compact object radius to unity: c = ro = 1 . For 
a maximally-spinning BH with a characteristic mass M = 3M 
dMacFadven & Wooslevl 1 1999b . our unit of length, in physical 



stellar envelope actively collimates the jet, we choose / such 
that it decreases with distance as / = (rlr m y y l 2 , giving the 
collimating shape ([1). Just outside the star, where pressure con- 
finement ceases, we allow the wall to start decollimating. To 
model the deconfinement we choose f{r) such that it reaches a 
minimum at r - r* and then, within a radial distance Ar = 0.5r», 
smoothly transitions to a logarithmic increase with increasing r. 
The exact manner in which the wall opens up is unimportant as 
we discuss in 50 

The code uses internal coordinates (xi, X2), which are uni- 
formly sampled with 1536 x 256 grid cells. These are mapped 
to the physical coordinates (r, 9) via r = 2.1 + exp(jq), and 
xj = sign( v F) \ y ¥\ 1 ^ 2 . In the lateral 9 direction, the computational 
domain extends from the polar axis X2 = 0, where we apply 
antisymmetric boundary conditions, to X2 — 1, where we place 
an impenetrable wall whose shape is determined by equation 
(0 with ¥ = I. In the radial direction, the computational do- 
main extends from r m = 3 (three radii of the compact staiQ) to 
'"out = 3xl0 n .Atr = r ou t we apply standard outflow boundary 
conditi ons, while at r — r m we apply jet injection boundary con- 
ditions (ITMN09I) . In detail, at r m we set the jet magnetic flux 
by specifying the radial magnetic field component according 
to (0, B I% Q = rr 2 , and set the mass flux ]mk into the jet by spec- 



ifying the plasma density po = 0.00132B 2 ., (s 



units, is GM/c 2 » 4.4 x 10 5 cm. 



2. Numerical Method and Problem Setup 



1.2 x 10 5 gcnr 3 
in cgs units for B r $ = 10 15 G) and poloidal 3-velocity vo = 0.5, 
both of which are independent of 9. 

We choose the rotational frequency of the central compact 
star to be D, = 0.25, which corresponds to a maximally spin- 
nin g BH, although we expect similar result s for other values of 



We have performed 
ativistic MHD code 



our simulations 
HARM 



ing recent improvements ( Mignone & McKinneyJ 12007 



using th e rel- ^ (IGammie et aU 120041; iMcKinnevl 12005b . We start the sim 



2003 



(Gammie et al. 

McKinnev & Gammiel l2004t IMcKinnevl l2006lh includ 



Tchekhovskov et al.ll2007t ItMN08UtMN09I) . We initialize the 
simulation with a purely poloidal field configuration in which 
the shape of each field line is given in polar coordinates by 



9(r) oc r 



-v/2 



0.5 < v < 1. 



(1) 



TMN08 showed that the above power-law shape corresponds to 
a power-law profile for the confining pressure: 



p oc r a , a — 2(2 — v). 



(2) 



We discuss our choice of v in {0 In the simulation, the jet 
shape (0 is maintained until a distance r = r», where r« rep- 
resents the radius of the confining stellar envelope. Beyond r», 
we allow the field to decollimate. 

In order to arrange for a smooth transition at r„, we suitably 
choose the flux function W(r, 6) of the initial magnetic field, 



ulation with a non-rotating star and smoothly turn on the rota- 
tion. This generates a set of outgoing waves as the initial purely 
poloidal magnetic field develops a helical structure in response 
to the rotation. In the steady-state solution, which quickly es- 
tablishes behind these waves and which we study below, energy 
flowing out of the compact star is the sum of (a) Poynting flux 
Jem, set by the rotation rate O and the field strength B r $, and 
(b) kinetic energy flux ;'ke, set by po and vq, Poynting flux is 
converted to kinetic energy as plasma accelerates along a field 
line and their ratio defines the local magnetization of the field 
line, cr = Jem/ 7ke- As the field line accelerates, magnetization 
cr decreases and Lorentz factor y increases in such a way that 
/i, the ratio of the total energy flux Jem + Jke to the mass flux 
7ma, is conserved, i.e., 



A* 



Jem + Jke 
7ma 



y(cr + 1) = constant along each field line, 



¥(r,0) = (1 - cos 6)1 f\r), 



(3) 



where we note that, by definition, T remains constant on field 
lines, and the total poloidal magnetic flux enclosed inside a po- 
lar angle 6 at radius r equals 27r v F(r, 9). The function f(r) essen- 
tially sets the opening angle of the wall and we adjust this func- 
tion so as to match our requirements. Inside the star, where the 



(4) 

where Jke = T7ma- Since cr > 0, the value of fi determines the 
maximum possible Lorentz factor of the field line. The Poynt- 
ing flux vanishes at the jet axis, so /j.{9 — 0) = yo = (1 - v^y 1 ^ 2 , 
and u inc reases away from the axis roughly proportional to *P 
(TMN09). reaching its maximum value, fig a 1150 (for our 
choice of boundary conditions), at the jet boundary. 



It is a matter of convenience where we locate the inner surface at which we 
inject the jet. We have verified that the solution is insensitive to the choice if 
the flow is sub-Alfvenic at this surface. 
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Figure 1 : Meridional jet cross-section, showing color-coded logarithm of the Lorentz factor y overlaid with poloidal field lines (thin solid lines corresponding to 
T 1 ' 2 = 0.1,0.2, . . . , 1). Thick solid lines show the position of the wall and dotted lines show the half-energy field line. [Panel a]: Jet confined continuously by a 
wall out to a large distance (model Moo). [Panel b]: Jet confined by a collimating wall until r = r, = 3 X 10 4 (model M4, where r, is shown by the horizontal dashed 
line), after which the wall opens up and the jet becomes deconfined. At this point, the jet separates from the wall and is surrounded effectively by vacuum. Once the 
jet becomes deconfined, y increases abruptly and the jet opening angle increases slightly. 



In the context of GRBs, we are interested in a jet which 
asymptotically has a Lorentz factor of order a few hundred. 
Also, we assume that the jet energy goes into accelerated parti- 
cles, and thereby to prompt GRB radiation, via an MHD shock. 
In order to have efficient conversion of jet energy into ther- 
mal energy at the shock, it is known tha t the jet fluid must 
achie ve <x < 1 just prior to the shock ( IKennel & Coroniti , 
19841) . These conditions require p ~ 10 3 and motivate our 
choice of po ~ 10~ 3 . By choosing this value of p, we design 
our simulation such that, if the jet achieves a Lorentz factor of 
a few hundred, it will have cr < 1 . 



3. Numerical Results 

We have run several simulations with the parameter v set 
to 3/4, and with different values of the confinement ra- 
dius r„ covering the likely range of progenitor star radii in 
GRBs. The choice y = 3/4 is motivated b y the wo rk of 
McKi nnev & Naravan d2007albl) and especially TMN08, who 
showed that this value corresponds to a confining pressure pro- 



file, p oc r~ 5 ^ 2 (eq. which is reasonable for a jet-shocked stel- 
lar envelope, e.g., see hydrodynamic simu lations of rel a tivisti c 
jets injected into GRB progenitor stars bv lZhang et al.1 (12003b . 
We have also run models with v = 2/3 (i.e., p oc r~ 8 ^ 3 ), and the 
results are not qualitatively different. 



As a baseline model, we first consider the case r„ — > oo, i.e., 
a model in which the jet is continuously collimated by a wall 
out to an arbitrarily large distance. We refer to this as model 
Moo. FigureQJ shows a meridional cut through the steady-state 
solution we obtain at the end of the simulation. In this model, 
the jet y increases steadily and continuously and becomes of 
order a few hundred at large distance, while the jet opening 
angle decreases continuously with increasing distance. 

Consider next a model with r* = 3 x 10 4 (= 1.3 x 10 10 cm 
for a 3M black hole). In this model, which we refer to as 
model M4, the wall collimates exactly as in model Moo until 
r ~ r„, but the wall then smoothly opens up (as described by 
eq.[3]l, allowing the jet to decollimate. The steady state solution 
is shown in Figure QJ). Just beyond r = r„ soon after the jet 
loses confinement, we see that (a) the Lorentz factor y increases 
abruptly by a factor ~ 10 and (b) the opening angle of the jet ; 
shows a slight increase and then "freezes out," hardly changing 
thereafter. 

Figure [2] shows the radial dependence of various quantities. 
In addition to the models Moo and M4 already discussed, we 
consider two additional models: model M3 with r* — 3 x 10 3 
(1.3 x 10 9 cm) and model M5 with r, = 3 x 10 5 (1 .3 x 10 1 1 cm). 
For each model, we show results corresponding to the "half- 
power" field line, i.e., the field line for which half the jet power 
is carried by field lines inside of this line and half by field lines 



3 



b 



le3 



le2 



10 



0.1 
1 

0.1 

le-2 



® 



10 r 



1 r 



0.1 V- 



le-4 fe- 
le-5 



•i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ti 




M3(r^=3e3) - _ 

, x M4(r =3e4) 

_W M5(r*=3e5) 

I I I I I II I I I I I I I I II 




M 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I'M 1 1 1 1 1 1 1 1 M 




3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I I 




~i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i 



0123456789 
log 10 r 



Figure 2: Dependence of various quantities along the half-energy field line for 
models Moo, M3, M4, and M5. These models differ only in the jet confinement 
radius, i.e., the radius of the star r*. [Panel a]: Lorentz factor y, magnetization 
o", and the total energy flux /i. Model Moo shows continuous smooth acceler- 
ation, whereas the deconfined models M3-M5 show a "burst" of acceleration 
soon after the jet loses confinement. Shown with color stripes is the analytic 
solution for a fully unconfined jet (eq.|5J. It reproduces our numerical solutions 
to within 20% (the width of the stripes). [Panel b]: Field line opening angle Oj. 
This decreases continuously while the jet is confined by the wall, but it hardly 
changes outside the star. [Panel c] : The product of the field line y and the full 
jet opening angle &j. In model Moo, we have y&j < 1 at all r. However, in 
models M3-M5, y& j increases abruptly soon after the jet becomes unconfined 
and reaches values ~ 10—30, as appropriate for GRB jets. [Panel d]: Magnetic 
pressure p m . This decreases smoothly with increasing r so long as the jet is 
confined (r < r»). However, once the wall opens up, p m drops abruptly. 



outside. 

Figure [2^ shows as a function of r the total energy flux [i 
(which is constant and the same for all models), the Lorentz 
factor y, and the magnetization parameter cr. While model Moo 
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Figure 3: Angular distribution of y (panel a), magnetic pressure (panel b), and 
power output per unit solid angle (panel c) in our fiducial model M4 as a func- 
tion of 6j, shown at various distances (see legend). Thick dots indicate the 
position of the half-energy field line. The collimation angle of this field line 
is nearly constant, 8j ss 0.024, once the jet becomes deconfined. However, its 
Lorentz factor increases dramatically. Crosses indicate the full jet opening an- 
gle &j fj(3}. Once the jet leaves the star, magnetic pressure and power output 
along the outermost field lines drop essentially to zero. Thus, outside the star, 
the jet is surrounded by a vacuum. 



shows continuous smooth acceleration out to large radius, each 
of the three deconfined models, M3, M4 and M5, exhibits an 
abrupt increase in y soon after the jet loses confinement, with 
most of the acceleration occurring between about r* and 2r*. 

Figure [2}} shows the dependence of the angle 9j between the 
poloidal component of the half-power field line and the jet axis. 
We use 9j for the opening angle of a given field line and ®j 
for the opening angle of the entire jet; for the latter we use the 
field line that contains 90% of the jet's power output. Figure^ 
shows that, so long as the jet is confined (r < 6j decreases 
with distance as r~ 3 / 8 (eq. Q]with v = 3/4). However, once 
the jet crosses r«, the opening angle of the field line freezes 
and remains nearly constant out to arbitrarily large distance. 
Therefore, it is the opening angle of the wall (eq.0 as evaluated 
at r — r* that determines the asymptotic value of the jet opening 
angle. We note that the opening angles of our jets are similar 
to those found in h ydrodynamic simulations of GRB jets (e.g., 
IZhang et aL . 2003 ). This is perhaps not surprising since we 
calibrate the shape of the confining wall in our models by the 
shape of the jet boundary in the hydrodynamic simulations. 
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After the burst of acceleration between r„ and 2r„, our de- 
confined jets subsequently accelerate very weakly. This is 
similar to unconfined magnetically-accelerated outflows which 
are kn own to accelerate only log arithmically at large dis- 
tance dBeskin et al.l Il998l ITMN09I). This is in contrast to 



thermally-accelerated outflows which continue to accelerate ef- 
ficiently until they reach their maximum Lorentz factor (e.g., 



Zhang et all 120031: iMorsonv et al.1 120071) . 



Once an initially confined jet is allowed to become decon- 
fined, its structure rearranges and begins to resemble a fully un- 
confined jet. Figure^ shows as colored stripes the ap proximate 
analytic solution for unconfined outflows obtained by lTMN09t 



r = 2 



y CIR 
-log 



y/3 


( \ 






> Jc = 




J 


K sin 2 8j , 



1/3 



(5) 



Here R = r sin Gj is the cylindrical radius, log is the natural 
logarithm, and Gj = (1 - v/2)(2 ( P) 1/2 (r H ,/r in )- v/2 is the angle of 
the half -power field line, evaluated at the stellar surface (equa- 
tion (0 with *P « 0.7 and / = (r,/ f"in) 2 ). The agreement 
of the analytical solution with our numerical results is better 
than 20%. This agreement suggests that, once deconfined, our 
jets rapidly relax to the equivalent unconfined monopole solu- 
tion for the same angle Gj by converting their excess magnetic 
energy into kinetic energy. There is, however, an important dif- 
ference. In a pure monopole solution, most of the jet energy 
flows out in the equatorial region, whereas in our jet models, 
because of the prior collimation, the energy flow is confined to 
within <8j < 

Figure |2j; shows the radial variation of the product y®j, 
where ; is the angle of the jet boundary (the field line with 
¥ = 0.95, which encloses 90% of the jet's total power). For 
the continuously collimated jet in model Moo (dotted line), y®j 
is of order unity throughout the acceleration zone, as noted by 
ko3 decreasing to a much smaller value at large distance. The 
unconfined jets in M3, M4 and M5, however, behave very dif- 
ferently. Once these jets become unconfined, y®j increases 
rapidly from =a 1 to ~ 10-30. An angular power-density- 
weighted average of yGj gives values within 30% of this esti- 
mate, which shows that the details of how we select y and ®j 
are unimportant. (A jet break is expected only if a large fraction 
of the power has y6 » 1 ; the power-density weighting accounts 
for this.) 

What determines the value of y®j for a jet? To find this 
out, let us rewrite equation (|5} in terms of magnetization, using 
equation ©: 



y®j = 2 3/2 (®j/ej)<T l/2 \og 



1/2 



7c 



15a- 



1/2 



(6) 



where in the last equality we assumed characteristic values for 
jet parameters in model M4 (r = 10 9 , 6j = 0.02, Q, = 0.25, 
fi = 10 3 , ®j/0j - 1.5). From this equation we clearly see that 
jets with subdominant magnetic fields at large distance, cr < 1, 
are limited to y®j < 15. 

We note that higher values of y®j correspond to deconfined 
models with smaller values of r„. From the simulation results 



shown in Figure[2j;, we obtain a scaling 

y®j <* r; - 22 . 



(7) 



We can use the analytic solution of fully unconfined jets to un- 
derstand this trend. As we show in the Appendix, for mod- 
els M3-M5, which are asymptotically roughly in equipartition 
(y ~ fi/2 at large distance), the analytic solution (|5]l gives: 



y — /u/2 : ®jOcr t v ^ 2 , y cc y ltL r v J q ', y®jOc/j. l/ 'r t 

(8) 

This implies y®j oc r~ - 19 for v - 3/4, in good agreement with 
the numerical scaling (0. In the limit of high magnetization, 
i.e., cr » 1 and y <K fi, we obtain instead (see the Appendix) 



l/2„v/4 



,l/2„-v/4 



y « fi : ®j oc r - v '\ y oc /'X , y®j oc fi l/s r; v,a . (9) 

Therefore, smaller values of stellar radius r» and larger val- 
ues of fi give larger values of y®j- However, in this limit, 
the asymptotic value of cr will be » 1. Such jets cannot ef- 
ficiently convert their Poynting energy flux into particle ther- 
mal energy via an MHD shock and therefore it is not clear if 
they will have sufficient energy in accelerated electrons to pro- 
duce the prompt y-ray emission observed in GRBs (see end of 
©. Other particle acceleration mecha nisms, e.g., magnetic re- 
connection or plasma instabilities (e . g . . |B eloborodovl I2009D . or 
different mechanism s, e.g., inverse Compton scattering of ther- 



mal photons (e. g.. iBrodericki 120051) , may be able to circumvent 
this limit. 

Finally, Figure [2jl shows the profile of magnetic pressure 
along the half-power field line. In the confined region of the 
flow, the pressure varies roughly as r~ 5 ^ 2 , leveling off to ~ r~ 2 at 
larger distance^ However, as the jet emerges from confinement, 
poloidal field lines rarefy laterally and the pressure decreases by 
over a factor of 10 (see models M3-M5). The drop in pressure 
creates a strong longitudinal pressure gradient which is the rea- 
son for the abrupt acceleration of the je t. This is an extreme 
example of the "magnetic nozzle" effect (IBegelman & Hi 1994 
TMN09I) in action. 

Figures [3^,b,c show for model M4 the transverse profiles of 
y, magnetic pressure p m , and power output per unit solid angle 
dPjdix) at various distances. We see that the abrupt acceleration 
and pressure drop experienced by the half-power field line (in- 
dicated by thick dots), described earlier in Figures |2^,d, is rep- 
resentative of most field lines. The outermost field lines near 
the edge of the jet (Gj > ®j) accelerate enormously during de- 
confin ement but they contain little power (< 10%) (ILvubarskv 
20091) . Along these field lines the pressure drops essentially to 
zero once the jet emerges from the star, so these outer field lines 
are surrounded by vacuum. 

We have included a fairly gentle opening up of the wall in 
the simulated models (e.g., see the shape of the wall in Fig.QJ). 
However, because the jet moves ultrarelativistically and all sig- 
nals are strongly beamed forward, even small changes in ge- 
ometry are equivalent to nearly complete deconfinement, e.g., 



2 The pressure profile in the slower magnetized sheath that surrounds the jet 
in a real system continues to follow the p oc r~ 5 / 2 dependence 1TMNO8). 
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the jet surface in Fig.QJi clearly separates from the wall and the 
pressure at the jet surface (Fig.[3j;) goes to zero. We have con- 
firmed that our results are insensitive to the precise shape of the 
opening-up of the wall. Since the wall is causally disconnected 
from the flow, a more strongly diverging wall shape would not 
affect the results. 

It should be noted that the numerical simulations described 
in this paper are extremely challenging since at high magnetiza- 
tion the equations of motion are very stiff and hard to solve. Ad- 
ditionally, it is a challenge to obtain accurate solutions that ex- 
tend over 10 orders of magnitude in distance and in which ideal 
MHD conserved quantities are preserved along field lines to 
better than 15% ( K09t TMN09I) : our numerical models achieve 
this accuracy and are well-converged. To ensure that we have 
sufficient angular resolution near the pole and the wall, and 
sufficient radial resolution to resolve abrupt radial changes of 
quantities where the wall opens up, we have rerun model M4 
with a locally 5 times higher angular resolution near the pole 
and near the wall, and also with a 5 times higher radial resolu- 
tion between 0.5r* and 4r». We obtained less than 2% difference 
in asymptotic values of y and 6j for most of the jet cross-section 
(0.2 < 9jl®j < 1). After this work was posted on the a rchives, 
our results were confirmed bv lKomissarov et al. (2010) using a 
completely independent code. 

4. Discussion 

GRB afterglow radiation is produced when a relativistic jet 
ejected from a collapsar interacts with an external medium and 
decelerates. The afterglow phase typically starts at the end of 
the prompt gamma-ray emission, say at a time to ~ 30 s in the 
observer frame. The Lorentz factor y and jet opening angle 
&j discussed in the previous sections refer to the jet proper- 
ties at this time. During the subsequent afterglow phase, the 
Lorentz factor y of the jet decreases with time. In the case of 
a uniform external medium and adiabatic evolution, y varie s 
with expansion distance as r 3 ^ 2 (iBlandford & McKee , Il976h . 
whic h corresponds to a variation with observed time as r 3 ^ 8 
(e.g., Sari et all 1998 : Piranl200"5 ). The opening angle of the jet 
&j, however, does not change. Therefore, during the afterglow 
phase of a GRB, the quantity y&j varies with observer time as 



y&j « (t/t y 3/s (y&j)o 



(10) 



where the subscript indicates values measured at the end 
of the prompt emission stage. When y@j falls below unity, 
there is an achromatic break in the afterglow light-curve and 



the observed flux falls more steeply with time (iRhoadsl 11999 
1999). This "jet break" has been seen in a num 
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ber of GRBs, typi cally abou t one to ten days af t er the initial 
prompt emission (IFrail et al. . l200lt ICenko et all l2010b . For 
these GRBs, according to equation (fTUl i. we require an initial 
(y&j)o ~ (1 day/30 s) 3/8 ~ 20 to explain the observations. 
This constraint has to be satisfied by any model of a GRB jet. 

We find that both confinement by the collapsar's stellar en- 
velope and free propagation outside the star are needed for an 
MHD jet to have (i) the required value of (y&j)o, (ii) the large 



power inferred in long-duration GRBs, and (iii) efficient con- 
version of electromagnetic to kinetic energy (cr < 1 at large dis - 
tance). Purely unconfined jets do give (y®j)o » 1 ( TMN09I) . 
However, they have cr » 1 over most of their volume, which 
is undesirable, and they have too little power within the rele- 
vant opening angle where cr < 1. Purely confined jets are much 
better at focusing the power of the central engine into a colli- 
mated beam with cr < 1, but they give (y®j)o ^ 1. Only jets 
that are first confined and then deconfined can explain the most 
energetic long-duration GRBs with achromatic jet breaks. Such 
jets abruptly accelerate as soon as they emerge from the stellar 
envelope and develop a typical value of (y&j)o ~ 20, while at 
the same time achieving cr < 1 (Fig. [2j». Remarkably, the stel- 
lar confinement radius we need to produce the observed value 
of (y&j)o is r„ ~ 10 4 ~ 5 or a physical radius ~ few x 10 10 cm, 
which agr ees well with the estima ted radii of GRB progenitor 
stars (e.g.. lHeger & Langer , 2000). 

Achromatic jet breaks have been observed in 20% of the 
Swift GRBs, but the actual number is much higher (~ 50%) 
since obser vations of most GRBs ceas ed before the expected jet 
break time (iKocevski & Buflerl 120081) . As we have shown, cur- 
rent observations are putting stringent constraints on GRB mod- 
els. If future observations by the Fermi Obse rvatory routinely 
find y > 10 3 dKumar & Barniol Duranl 120091) and (&j) Q > 0.1 
dCenko et all 120101) for many GRBs, we would need to explain 
why these jets have {y®j)o > 100. The numerical factor in 
equation © is a logarithmic factor and is unlikely to be larger 
than 10-20. Therefore, either we must accept that cr » 1, 
i.e., GRB jets are Poynting-dominated and somehow manage to 
convert a large fraction of their energy to prompt gamma-rays, 
or that MHD is not the appropriate framework for understand- 
ing GRB jets. 

While we have focused primarily on GRB jets, our results 
also have some bearing on relativistic jets in XRBs and AGN. 
In order to produce collimated and efficient (<x < 1) jets, the 
accretion disks in these systems must collimate the outflow- 
ing Poynting flux from the central BH. Strong winds are natu- 
rally present in advection-dominated accretion flows (AD AFs , 
Naravan & yJ 1 1994 1 19951: iNaravan & McClintockl Eooil). so 



one expects to see well-developed jets whenever XRBs and 
AGN are in the ADAF state. The exact properties of the emerg- 
ing jet will presumably depend on the distance to which effi- 
cient wind-collimation operates. 

We note that our simulations assume a fixed shape for the 
confining wall. However, as the GRB progresses and the 
jet disrupts the star, the opening angle of the jet is likely 
to increase. Our model does not take this effect into ac- 
count. Also, our assumption of a rigid wall as the jet con- 
fining agent eliminates instabilities and mixing at the jet in- 
terface. These effects should be taken into account in future 
work, although the problem of a magnetized ultra-relativistic jet 
interacting with a surrounding progenitor star envelope is nu- 
merically extremely challenging and remains unsolved. Future 
work should also stud y jet stability to non-axisymmetric modes 
(INaravan et all 2009), general rela t ivistic effects of the central 
spinning BH dTchekhovskov et all |2010b . how GRB engines 
become threaded by the ordered magnetic fields needed for 
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MHD jets dMcKinnevl l2006t iMcKinnev & Blandfordl 120091) . 
and should model radiation processes in order to make de- 
tailed comparisons with observed afterglow light curves and jet 
breaks. 
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Appendix A. Analytic scaling of yQj in jets 

In this section, we analytically derive relations for y®j. For 
this, we substitute &j oc r» v ^ 2 (see eq. Q} into equation (0 and 
neglect the logarithm: 



fi-y 



(A.l) 



Now, we take a logarithm of both sides of equation ( 1A.1I ) and 
take a differential: 



^dy dfj. - dy dr 
y fJ--y r* 



(A.2) 



Our models M3-M5 are asymptotically roughly in equiparti- 
tion: y =s fi/2 at large distance. Therefore, we substitute 
p - y = y = fx/2 in the result, 



4^Z _2— = v— 

y fi r. ' 



(A.3) 



Integration of this equation gives: 



y^H/2: ; -ocr; v/2 , y oc /j. 1/2 r Y J 4 , y&j oc // 1/2 r; v/4 - 

(A.4) 

This implies y®j oc r~ 019 for v = 3/4, in good agreement with 
the numerical scaling (0. In the limit of high magnetization, 
i.e., <x » 1 and y <k fi, we approximate ji - y as fi in equation 
( IA.U and obtain instead 

y«,u: ®jocr; v/2 , ycc^^r v J\ y®j oc n l/i r ; v/6 . 

(A.5) 



7 



